Entanglement and decoherence in spin gases 
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We study the dynamics of entanglement in spin gases. A spin gas consists of a (large) number of interacting 
particles whose random motion is described classically while their internal degrees of freedom are described 
quantum-mechanically. We determine the entanglement that occurs naturally in such systems for specific types 
of quantum interactions. At the same time, these systems provide microscopic models for non-Markovian 
decoherence: the interaction of a group of particles with other particles belonging to a background gas are 
treated exactly, and differences between collective and non-collective decoherence processes are studied. We 
give quantitative results for the Boltzmann gas and also for a lattice gas, which could be realized by neutral 
atoms hopping in an optical lattice. These models can be simulated efficiently for systems of mesoscopic sizes 
(N ~ 10 5 ). 
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We study the entanglement properties of spin gases. A spin 
gas is a system of interacting spins (or qubits) where the cou- 
pling strengths between the spins are stochastic functions of 
time. A system that could serve as a textbook example of a 
spin gas is the semi-quantal Boltzmann gas, where each par- 
ticle carries an internal (two-level) quantum degree of free- 
dom. During a collision of two particles, the internal degrees 
of freedom interact and can become entangled. The statis- 
tics of the collisions, described by kinetic gas theory, leads 
to randomly fluctuating coupling strengths between the spins. 
A similar situation arises in systems where the gas particles 
do not move freely in space, but are confined to lattice-sites, 
between which they can classically hop with a certain prob- 
ability. It is an intriguing question, how the evolution of the 
quantum state of the system is determined by the underlying 
classical thermodynamics of the gas. What kind of entan- 
glement is created in the gas and at which rate? How is the 
equilibrium state of the gas characterized in terms of its en- 
tanglement? In this paper we will give an answer to these and 
to other questions, which make the study of spin gases inter- 
esting both from the perspective of thermodynamics and of 
quantum information. 

Spin gases differ from spin lattices in that the coupling 
strengths have no translational symmetry and evolve in time. 
Spin gases are more closely related to spin glasses, which have 
random, albeit static, couplings between the spins. Although 
there has been much recent work investigatinge.g. the role of 
entanglement in quantum critical phenomena HHSSB IB 0] > 
there is little theoretical work studying disordered quantum 
systems such as spin glasses and spin gases from a similar 
perspective (see, however, 01). 

We give a simple, and yet realistic, collision model for the 
particles carrying the spins, and study the quantum mechani- 
cal states that emerge from the dynamics of such a spin gas, 
which we also refer to as a semi-quantal gas. In full gener- 
ality, this problem seems intractable for various reasons: the 
description of a many -body quantum state usually requires ex- 
ponentially large resources; strong interactions do not allow 



for a perturbative treatment; random interactions prevent the 
appearance of symmetries and the corresponding reduction of 
the effective number of degrees of freedom; finally, the re- 
striction to low-energy eigenspaces, suitable for the study of 
ground-state or low-temperature properties, can not be ap- 
plied here to the study of dynamics. State of the art numeri- 
cal methods — such as density matrix renormalization group 
— are limited to systems of moderate size (up to a few hun- 
dred particles) with bounded amount of entanglement Gl (e.g. 
non-critical spin chains) or an amount of entanglement that 
scales at most with the surface of the block of spins 0] (e.g. 
some 2-dimensional spin lattices). In disordered quantum sys- 
tems with random interactions, like the spin gas studied here, 
entanglement will typically increase with the volume of the 
block, which implies that these numerical methods can not be 
applied. Nevertheless, for Ising (or, more generally, commut- 
ing two-body) interactions we can compute the full dynam- 
ics of the many-body system exactly and efficiently. Semi- 
quantal gases are not only toy models of theoretical interest, 
but could be experimentally realized even with present-day 
technology and existing setups. 



In this paper we characterize the states that arise in spin 
gases, and calculate their expected entanglement properties. 
We point out that we describe the system by a pure state, albeit 
with random coefficients, and calculate expected properties 
by averaging them over different realizations of the state. The 
formalism covers a range of different models of spin gases and 
allows us to identify many of their quantum features. For in- 
stance, we can describe the quantum dynamics of a particular 
set of (probe) particles prepared in various initial states. The 
gas provides a microscopic model for non-Markovian deco- 
herence 11(11 that can be treated exactly even for large system 
sizes. The paper is structured as follows: we first introduce 
the formal framework to study spin-gases, describe general 
entanglement properties, and then treat two specific models 
of a spin gas. We present analytical results for semi-quantal 
Boltzmann gas, and numerical results for a lattice gas. 
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Formal preliminaries. We consider a situation where N 
particles move along some classical trajectories r k (t), while 
their quantum degrees of freedom interact according to a 
distance- and time-dependent Hamiltonian 

k<i 

where g(rk(t),ri(t)) is some function depending on the par- 
ticular two-body interaction. We restrict ourselves to specific 
types of interactions, namely those where all ij( fci ) commute. 
Consequently, we find that after a time t the initial state |\&o) 
evolves to 

\*t) = U t \*o} = l[uW(vM(t))\*o}, (2) 

k>l 

with U^ kl \ip k i(t)) = e -Wt)ff ( *° and 

m(t) = I g(r k (t'), ri (t'))dU '. (3) 
Jo 

At a time t, the quantum state is fully determined by the 
N(N — l)/2 phases (fiki(t), which in turn are determined by 
the interaction history of the N particles. Each phase can be 
interpreted as a matrix element T k i = ip^i of an adjacency 
matrix T(t) defining a weighted graph. Thus, many entangle- 
ment properties of the state can be expressed in simple graph- 
theoretical terms. 

We now focus our attention on the case of an Ising-type 
interaction with = |11} H (11| 0. We further as- 

sume that all particles are initially prepared in the internal 
state |+) = 1/V2(|0) + |1)), and interact only when they 
collide, i.e., g(r k i) = 9o forr M (i) = \r k (t) - ri(t)\ < d Q and 
g( r ki) = otherwise. These restrictions are made for sim- 
plicity, but similar methods can be applied to general commut- 
ing H^ kl \ arbitrary g(r k (t),ri(t)) and pure separable initial 
states (without additional overhead). 

The evolution of the initial state \^o) = |+)' 8Ar 
can be straightforwardly described in the standard basis 
{|0),|1>}®", 

U t \+f N = 2-f J2 U * I s ) = 2 ~ f E e^ s ' r(t) ' s |s) , (4) 

s s 

where the sum is carried out over all iV-digit binary vectors 
s, i.e., over all 2 N different combinations of zeros and ones. 
We note here that all the time dependence is in the adjacency 
matrix T(t) of the graph. The parametrization of the quantum 
state in terms of a weighted graph that summarizes the "col- 
lisional history" of the gas is both intuitive and useful for our 
computations. 

Many properties of this global pure state can be understood 
in terms of the reduced density matrices of its subsystems. 
Since the unitary operations in commute with each other, 
the evolution of a set A of Na particles can be separated into 
two contributions. The first entangles particles within A and 
is determined by the block Taa of the adjacency matrix. The 



second contribution couples the subsystem A to the rest B of 
the system through the off-diagonal block Tab- The effect of 
the latter can be obtained by tracing out the set of particles B 
from the state |^ t ): 

2 N -l 

PA = ^tr B £ e ^(s.r.s-s'.r.s') ^ 

= air E(2^E el(SA ^ ) ' r ^- SB )l s ^l(5) 

The second equality is obtained by writing \s) = \sa) \sb), 
and the tilde in pA indicates that interactions within subsys- 
tem A are not taken into account {Taa is set to zero). Clearly, 
the values of the block Tbb do not affect the properties of 
either pa or p~A ■ In the standard basis, each off-diagonal ele- 
ment ("coherence") of the initial state is decreased by a factor, 
Ps A s' A (t) = C SAS ' A p SAS ' A (0), while diagonal elements remain 
untouched. The multiplying factor can be conveniently writ- 
ten as 

N B 

c sasW = e '^' (s - s ^ Tk J] cos ^ - s 'a) ■ r k ] (6) 

fe=l 

where we have defined the Na -dimensional vector (Tk)j = 
T k j for each particle k 6 B. In this form, we see that the 
total effect of the interactions with particles in B on a partic- 
ular coherence C SaS ' a of pa can be obtained by multiplying 
the effects of each individual particle in B. More succinctly, 
if p^ is the state of the subsystem due to the sole effect of 
particle k E B, then the state pA is obtained (up to normal- 
ization) by the Hadamard product of all {Pa fk=i wr i tten m 
the standard basis, that is by their component- wise multiplica- 
tion. This observation can also be understood within the con- 
text of Valence Bond Solids (VBS) as was recently shown in 
The decomposition into Hadamard products allows one 
to read off the matrix elements of pa from the adjacency ma- 
trix, but most importantly it signifies that one can efficiently 
compute reduced density operators pa of small subsystems 
A, even when the size N of the total system is essentially ar- 
bitrarily large. The computational effort scales only linearly 
with \B\ in contrast to the general case where the computa- 
tional resources to calculate pa scale exponentially with \B\ 
(because the partial trace has to be performed over all 2' s ' 
basis states in B). 

The time dependence of the quantum state of the system 
in terms of T(t) @, together with the efficient method (|6jl to 
compute the state of (small) sub-systems are crucial proper- 
ties that allow us to study spin gases. For a complete char- 
acterization of the dynamics, one still needs to find, for each 
particular gas model, the behavior of the stochastic function 
r(i). This task amounts to assigning a probability Pr(t) to 
every collisional history T(t). However, in order to calculate 
the evolution of average properties, it is enough to have the 
probability distribution of the T's at every given time pt(T). 
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Entanglement and decoherence. While entanglement 
of bipartite systems is rather well understood, entanglement 
properties of multipartite systems are in general difficult to de- 
termine. However, for pure global states, we can get a broad 
picture of the entanglement in a multipartite system by con- 
sidering all possible splits of the set of parties in two groups 
(A and B), and analyzing their bipartite entanglement. For 
each of the 2 N ~ l bipartitions, the entanglement properties can 
be completely determined from the Schmidt decomposition, 
Y&t) = Sfc \f^k\k) A\k) b, where the Afc are the eigenval- 
ues of the reduced density matrix pa = tiB\^t}{&t\- The 
entropy of entanglement Sa = tr(p^ log 2 pa), i.e., the von 
Neumann entropy of the reduced density operator pa, pro- 
vides a suitable measure of bipartite entanglement, quantify- 
ing bipartite aspects of multipartite entanglement. However, 
as mentioned before, the calculation of reduced density matri- 
ces is, in general, very difficult if not impossible (exponential 
scaling in both Na and iVg). Even the simple criterion to as- 
certain whether the system is entangled or not (rank(/?A) > 1 
or = 1) may be impossible to check. For states |^t) which 
occur in the spin gases under consideration, many of these re- 
strictions do not apply. First, we can use (|6) to determine in 
an efficient way the density matrix pa of small subsystems A 
and hence calculate the entropy of entanglement with respect 
to all such bipartitions of the system. Interactions within A 
do not change the entanglement with the rest of the system 
and can therefore be ignored. We can also neglect the phase 
factors in C ss ' ©, since they can be cancelled by applying lo- 
cal unitaries V^j = Ski cxp(-i|ij m rj m ) on every qubit 
j 6 A. Second, we have a simple criterion for the presence 
of entanglement 12^1 : The state \^>t) is entangled with re- 
spect to the partition A — B iff the two groups are connected 
(i.e. an interaction between some particle in A and some par- 
ticle in B has taken place). We stress that since we can com- 
pute pa of small-sized subsystem, we can also compute quan- 
tities such as the multipartite Meyer- Wallach pure-state en- 
tanglement measure 11311 . which only depends on single -body 
density matrices, or the "correlation strengths" fl4ll for finite 
blocks, which depend on the reduced density matrices of the 
block and all of its reductions, or classical correlation func- 
tions like those used in generalized n-party Bell-inequalities 
Q. 

A different approach to describe global aspects of the mul- 
tipartite entanglement present in the gas is to ask whether or 
not entanglement can be created, or localized, between two 
arbitrary subsets of particles A1 — A2 by performing local op- 
erations on the other particles. Again, we find that for the class 
of states that emerge in the spin gas there is a very simple cri- 
terion to answer this question: entanglement between A1—A2 
can be created by local operations iff there exists an path be- 
tween A\ and A2 in the corresponding graph. The necessity 
of the condition is obvious, while the sufficiency follows from 
these two facts: (i) Collisions where particle k is involved can 
be undone (up to local operations) by measuring k in the z- 
basis {|0) , |1)}. (ii) If some particle i is connected to a par- 
ticle j via an intermediate particle k with generic phases ipik 



and ipkj, then a projection onto a state \x) k ^ {|0) , |1)} ef- 
fectively creates an entangling operation on the particles i and 
j (unless both are in one of the states of the standard basis). 
One can hence perform measurements (e.g. along the x-axis) 
on the particles found in the connecting path, z-measurements 
on the rest, and ends up with an entangled pair. Determining 
whether two particles in a graph are connected is known as the 
reachability problem, which can be solved by an algorithm of 
0(N 2 ). One can try to optimize this procedure to obtain the 
maximum average entanglement, called localizable entangle- 
ment. This quantity, which has recently been introduced in 
0], is bounded from below by the two-point classical corre- 
lation function, and from above by the entanglement of assis- 
tance. Both bounds can be readily calculated from the reduced 
density matrix pa x a 2 of the subsystems and can therefore be 
computed efficiently for the spin gas. 

Let us now adopt a slightly different point of view that will 
prove useful in explaining some central features of our mod- 
els. We focus our attention on a particular set of (probe) parti- 
cles, and study the map £ t (p) induced by their collisions with 
the other particles, which we call the background gas. On the 
one hand, the action of this map on the state &)j E A gives 
Pa, which provides us with the necessary information to cal- 
culate the entanglement between A and the rest of the system. 
On the other hand, the action of the map can be understood 
as a decoherence mechanism acting on general states of probe 
particles. 

A map is fully specified by its action on an operator-basis. 
In particular, it is straightforward to show that £ t (\s)(s'\) = 
C SS ' \s)(s'\, with C ss i given by 0, as long as the state ps of 
the background particles fulfills tr(|l) fe (l| ps) = 1/2 for all 
k £ B 1 27]. This means that in a possible physical realization 
of such decoherence studies one can relax the experimental 
conditions that would be necessary to keep the background 
gas in a coherent state. In addition, the Hadamard product 
structure of the maps has the same implications that we saw 
for states l28ll . Semi-quantal gases provide therefore an inter- 
esting model for decoherence in a mesoscopic environment. 
These models can be treated exactly and can be highly non- 
Markovian (i.e., show memory effects) in some regimes. 

From we find that if a part Al of the probe subsystem A 
does not have any common collisional partners with the rest 
A2 of the subsystem then C SA ^ A = C SAltSA1 C SA2 y A2 where 
sa = (sai,sa2)- Accordingly, we will say that the maps or 
channels acting on each part of the sub-system are indepen- 
dent (or uncorrected) 12911 . On the other hand, if two parts 
of the subsystem share collisional partners then some coher- 
ences will be nearly unaffected by such correlated collisions, 
while others will suffer an increased decay (as compared to 
uncorrelated collisions). For instance, consider the case were 
two probe particles (1 and 2) have a very similar collision pat- 
tern, i.e., T2j = Tij +Sj for all j, then it follows that co- 
herences associated with |01)(10| will only decay by a factor 
2~ Nb J2 Sb e lS ' SB , while |11)(00| will be "super-damped" by 

o-JVb i(2r 1 +<5)-s B 
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The above conclusions follow directly from the values of 
the adjacency matrix at a given time. We now proceed to study 
the explicit time dependence of the map, where we concen- 
trate on one central aspect: Markovian versus non-Markovian 
dynamics. Both effects can already be seen in the evolution 
of a single probe particle, and essentially correspond to the 
following two types of collision patterns that occur during a 
short time interval At: (i) Markovian: At every time step St 
the probe particle collides with a different particle and accu- 
mulates a small interaction phase Sip. The state of the probe 
particle will then decohere exponentially fast with the number 
of time steps k = At/St: \p 01 \ = [cos(<5<^/2)] fc = e~ At / T ' 
with r e w 85t/5ip 2 . (ii) Non-Markovian (or coherent cou- 
pling): In the time interval At a given gas particle has col- 
lided k times with the probe. The coherent addition of the 
interaction phase leads to a Gaussian type of decay: \poi\ = 
cos(fc<V/2) w e " A * 2/(2T 9 2) with r g = 2St/Stp. Addition- 
ally, if we assume no control over quantum or classical de- 
grees of freedom of the background gas, as in our the deco- 
herence analysis, we should average the effect of the maps 
£t(p) over all possible collision patterns at the given time: 
C s>a '(t) = J dT Pt (T)C s ^{T) where p t (T) is the probability 
that at a time t the adjacency matrix is T. So, the phase factor 
in (|5} suffers two types of averaging effects, one over collision 
histories, and the other one over all possible combinations of 
excited (i.e. state |1)) collisional partners, i.e. over all Sb- A 
complete characterization of the dynamics of the probe can be 
given in terms of the overall probability distribution p ZA (e). 
Here, e = za ■ Tab ■ sb is the phase in (0 for a given co- 
herence specified by binary vector za — sa — s^. One then 
finds that the coherence falls off essentially with the width 
(Tr(t) of the distribution as cxp(— ar(t) 2 /2) and acquires a 
phase $ ~ (e) . The two extreme regimes mentioned before 
correspond to or(t) oc s/t (Markovian) and or(t) oc t (Non- 
Markovian). We point out, however, that the dynamics will 
typically be very rich showing combinations of both effects 
and further non-trivial features such as finite-size effects. 

Boltzmann gas. We consider a dilute ideal gas of N parti- 
cles in thermal equilibrium with a mean free-path comparable 
to the size of the enclosing volume. The statistical state of 
the gas is fully specified by the density n, the volume V, and 
the temperature T. We assume Stosszahlansatz (or molecu- 
lar chaos) and hence take a homogeneous and uncorrected 
spatial distribution of the particles (density), and an uncorre- 
cted velocity distribution. The latter is given by the Maxwell- 
Boltzmann distribution and is characterized by the single pa- 
rameter a = y/ ksT jm, where m is the mass of the particles 
and ks the Boltzmann constant. We further assume a hard- 
sphere model for collisions between particles of diameter d 
and that at every collision particles acquire a phase inversely 
proportional to their relative velocity, tp kl = 'y/vki- We study 
the entanglement that arises in the system if at a given time 
t = the internal state of all particles is initialized to | +) . 

One could compute entanglement properties by direct sim- 
ulation of the Boltzmann gas. Here, however, we will focus 
on regimes where analytical results can be obtained, namely 



for large collisional phases and arbitrary times, or for arbi- 
trary phases and in the limits of short and infinite times. In 
what follows we will use the von Neumann entropy as well 
as the Reyni entropy as a measure of pure state entanglement. 
Renyi entropies, S q (p) = (1 — q)^ 1 log 2 (trp 9 ), for q > 
are known to be non-increasing functions of their parameter 
q, i.e. S q < S q i for q' > q. Moreover, in the limit q — ► 1 the 
Renyi entropy coincides with the von Neumann entropy, and 
therefore S q= 2 = — log 2 (trp 2 ) provides a lower bound to the 
von Neumann entropy S(p) = S q ^i > S g -2- 

In a regime of large collisional phases (j> ~ ja^ 1 ^> 1 
(i.e., large interaction constant or low temperatures ) we can 
assign to each collision event a random phase in [0, 2ir}. This 
already allows us to find the expected entropy for short times 
rt < 1, where r = nd 2 n(v r ) is the collision rate and (v r ) = 
y/WKT^nr)^ 1 is the mean relative velocity. For these 
short times a particle will typically collide at most once (with 
probability rt), and an expected entropy of J= J S(4>)d<p = 
2 - log 2 e. Hence, we find (Si) = rt{2 - log 2 e) + 0{{rt) 2 ). 
The entropy of entanglement between a block A of size Na 
and the rest of the system B can be readily obtained by count- 
ing the collisions that typically occur between particles of A 
andB: 

(Sa) w ^prf r<(2 - lo 8'2 e) for rt < 1. (7) 

For arbitrary times, we can use the convexity of the loga- 
rithm function to obtain a lower bound for the von Neumann 
and Renyi entropies: (Sa) > S q=2 {pA) = -(log 2 (tr/j^)) > 
-log 2 «tr&» = log 2 (E s ,,^(|C s ,xJ 2 »- From ©we 
notice that the coherence C SA>s i A only depends on the differ- 
encez^ = s^— s^. And in particular, the average (|C SA;S ' A | 2 ) 
depends only on the number Za of non-zero entries of za- 
Each particle k in B will contribute with a factor 1/2 to 
the product in © if it has collided with at least one parti- 
cle of the subset of A where za has non-zero entries, while 
a factor one appears otherwise. Since the probability that 
no such collision occurs is p ZA = exp(— HZa/N), on av- 
erage each term in the product will contribute with a factor 
1/2(1 - p ZA ) + p ZA = (1 + cM-rtZ A /N))/2, and taking 
into account combinatoric factors we arrive at 

mt)) > -io g2 E (^)d + e - rt ^) . 

(8) 

Numerical results for < 8 and arbitrary system sizes 
show that this lower bound is also a good estimate and de- 
scribes well the behavior of the entropic entanglement. This 
lower bound can be supplemented by the upper bound (Sa) < 
Na(Si). The latter follows from the sub-additivity property 
of the von Neumann entropy. 

In the short- or long-time limits we can simplify the above 
expression for the entanglement between two arbitrary parts 
of the Boltzmann gas. For short times HNa/ (A — 1) < 1, 
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we find 

NaNr 

(5 A ( i ))>-log 2 (l- i ^ yrt ) 



N A N E 



41n2(iV- 1) 



rt. 



which is consistent with the exact result @ for short times. 

In the long-time limit all particles will have collided with 
all other particles many times and accumulated phases ipki ^> 
1, independently of the collisional phase per collision. We re- 
fer to such a state as the equilibrium state. As times increases, 
rt 3> N, every term with za > in the sum approaches 
one (disregarding the binomial factor) while the Za = term 
is equal to two, 



(Sa) > 



log 2 



1 



1 

2^ 



N A N E 



The equilibrium state (rt — > oo) has thus the interesting fea- 
ture that | "too) is maximally entangled with respect to all pos- 
sible bipartitions, i.e. Sa ~ Na, provided that the total num- 
ber of particles TV in the gas is sufficiently large. This is a non- 
trivial statement especially in the case where both partitions 
are similarly large. We find that, for whatever bipartition one 
takes, the expected entropy of entanglement is at most a single 
bit away from its maximal value: Na > (Sa) > Na — 1. This 
result is in agreement with the findings of Page 11711 and sub- 
sequent work studying average or typical entanglement prop- 
erties of the whole set of multipartite pure states. Another 
remarkable property of the equilibrium state is that the local- 
izable entanglement between any pair of atoms approaches its 
maximum value of one e-bit as N increases. The proof is sim- 
ple. For any pair of atoms {i,j} with an interaction phase ifij 
we can find a third particle k with which both atoms have ac- 
cumulated a phase tp ik m ipj k m n. Then, we measure the 
remaining gas particles in the z -basis (to effectively decouple 
them) and particle k in the a>basis. A straightforward calcula- 
tion shows that the resulting state of particles {i, j} is always 
maximally entangled (20], no matter what the initial tp^ was. 
Whereas the dynamics of the entanglement strongly depends 
on the particular type and regime of the semi-quantal gas, the 
equilibrium state will, in general, be of the type we just de- 
scribed. In particular, it describes the state of the semi-quantal 
Boltzmann gas also in the regime of small collisional phases, 
and the lattice gas studied below. 

To conclude our analysis of the semi-quantal Boltzmann 
gas, let us give the short-time entanglement in the regime 
of small collisional phases. In this regime, we have to take 
into account the relative velocity of the colliding partners. 
A simple kinetic theory calculation shows that at short times 
(rt <C 1) the expected value for the squared modulus of the 
coherence is given by, 



tmrd 2 / dvp r (v)v (1 - \C V \ 2 ) = 1 - at, 



<|On 



where p r (v) is the relative velocity distribution and |C U | 2 = 
cos 2 (0(w)/2) is the squared modulus of the velocity depen- 
dent coherence. For small collisional phases -fa^ 1 < 1 the 



proportionality factor a can be approximated by 

1 3 f v 2 . 7 

a = 4-K 2 d 2 n(- / dv v 3 e~^ sin 2 (-M « 

v Ana 2 J y 2v 
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Following the previous reasoning we find 

/Q NaNb i fi a + \ a NaNb 

{Sa ^ > - AF^T l0g2(1 - 2 <} w 2hT2 



t. (9) 



we obtain the rate 



Thus, by writing a 
of entanglement generation in terms of the thermodynamical 
variables. Whereas in the low temperature regime entangle- 
ment grows with the rate of collisions (oc VT), here we find 
that entanglement generation is governed by the slow colli- 
sion events (larger phases) leading to the opposite behavior 
(a oc VT^ 1 ). 

Lattice gas. The quantum properties of the system are di- 
rectly linked to the classical statistical properties of the gas 
through r(t).Thus, in general, it is necessary to know the 
classical ?i-body distributions to give a complete description 
of the quantum state. For some gas models and regimes (like 
the Boltzmann gas presented above) correlations play a minor 
role and one can find analytical descriptions for single-particle 
phase-space distributions. In the remainder of this paper we 
study a lattice gas model that exhibits strong correlations. We 
use numerical simulations to recreate such classical correla- 
tions — or, in other words, produce different random real- 
izations of r(t). The model has the additional feature that it 
can be possibly implemented in a quantum optical system. It 
has already been shown that an optical lattice can be used to 
store ultra-cold atomic gases. The degree of control in these 
experiments is extraordinary. Various system parameters can 
be adjusted — from the structure of the underlying periodic 
lattice to the interactions between the atoms — , opening the 
door to a wide range of experiments and theoretical proposals 
ll6l fl8ifl9[EolE4iE3.E3ll . In particular, one can choose a pa- 
rameter re gim e where each lattice site is occupied by at most 
one atom Il8l EUl . The internal quantum state of the atom 
(e.g. two meta-stable hyperfine states) can be stored in coher- 
ent superpositions over long time-scales (few minutes), while 
coherent inter-atomic interactions have also been achieved by 
cold collisions fiUl . These interactions correspond precisely 
to the Ising-type chosen here. One can also find schemes 12411 
to induce a random (incoherent) hopping of the atoms from 
one site to its neighboring sites. Hence, we consider an MxM 
lattice containing N particles that randomly hop from site to 
site with a hopping rate rj, and have nearest-neighbor interac- 
tions with coupling constant g Q . 

We use the above lattice gas model to investigate the typi- 
cal multipartite entangled states that arise. An interesting ob- 
served effect is the creation of clusters, that is, sets of particles 
that are connected in the graph, such that entanglement can be 
localized between any of the constituent particles. One finds 
that the average number Nc of particles that form a cluster 
grows exponentially with time. Similarly, the expected value 
for the maximal distance i max between (localizable) entan- 
gled particles increases. After a finite time to, which only 
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depends on the filling factor v = NM~ 2 and hopping rate 
77, almost all particles are connected. That is, after this time 
to, all particles are pairwise (localizable) entangled and hence 
f max ~ m, suggesting a (bond) percolation phenomenon in 
the underlying graph. 




ni lit 

FIG. 1: Lattice gas with M = 20, g = 0.8rj. (a) Average entropy 
of entanglement for block sizes =4 (dark) and \A\ = 8 (light) 
and filling factors v — 0.9 (dashed) and v = 0.25 (solid), (b) Aver- 
age entropy of entanglement for two probe particles that are initially 
located at a relative distance t = (thin-dark) and £ — 8 (thick-light) 
and have hopping rates rj p = (dashed) and r\ v = \r\ (solid). 

Concerning block-wise entanglement, Figure ^a) shows 
the expected entanglement between a set of randomly cho- 
sen particles and the rest of the system as a function of 
time. Since all the density matrix elements are initially equal, 
|\T/o) = |+} 8 ' JV , and all coherences decay to zero with time, 
we find that the entanglement saturates to its maximum value. 
At high filling factors the entropy of entanglement undergoes 
periodic decreases every time the typical accumulated phase 
4>ik(t) with a nearest neighbor reaches 2ir, i.e., at t w k2n / 1 g 
with integer k. At lower filling factors we can observe that 
the entanglement is not additive when increasing the block- 
size, which indicates the presence of correlated collisions. To 
highlight this effect, we consider the case where two probe 
particles move freely in the lattice — with a hopping rate r\ v 
possibly different from the hopping rate of the background 
particles — and only interact with the gas particles that are 
in the same position El. FigOb) shows the entanglement 
of the probes with the rest of the system for different initial 
relative distances I between the two probes, and for differ- 
ent probe hopping rates r) p . Small probe hopping rates favor 
repeated collisions with the same particle — as opposed to 
independent collisions with different particles. As we saw 
before, this amounts to coherent addition of the interaction 
phase, which leads to a much higher rate of entanglement gen- 
eration. At I — and fixed probes (r} p = 0), all the collisions 
with gas particles are exactly correlated. From the above dis- 
cussion we know that under these circumstances the coher- 
ences of the type poi.io remain untouched while all others 
decay to zero, leading to a maximal entropy of S(oo) = 1.5. 
When the probe particles are given a non-zero hopping rate, 
they start suffering independent collision events which even- 
tually turn the reduced density matrix into a completely mixed 
one (S = 2). 

In the previous plots we assumed that the probe particles 
were initially prepared in the same quantum state as the back- 



ground particles. In the following, we study the decoherence 
of different entangled initial states of the probe particles due to 
the action of the background gas. In other words, we analyze 
properties of the average state of the probes, as opposed to av- 
erage properties of the resulting states. We consider the situa- 
tion where the probe particles are displaced at a constant speed 
v relative to the gas. We vary the speed v of the probes as well 
as the distance d between the probe particles. By doing so, we 
can distinguish two different effects, analogous to those stud- 
ied in the context of entanglement generation: (i) By decreas- 
ing the probe speed, we analyze the effect of multiple interac- 
tions with the same particle in contrast to interactions with dif- 
ferent (independent) particles, (ii) By increasing the distance, 
we turn from correlated to independent interactions of the 
probes with particles in the remaining system. Fig. |2ja) shows 
the decay of the entanglement in the Bell state \<f> + ) as mea- 
sured by the concurrence 1 25] in two extreme scenarios: (i) 
The probe particles are fixed in their initial positions (v = 0). 
(ii) The probe speed is chosen large enough (v /a ^> 77, where 
a is the inter-site spacing) so that a fixed value <p = 0.1 can 
be assigned to the collisional phase every time a probe par- 
ticle crosses an occupied site. These two scenarios illustrate 
the difference between Markovian and non-Markovian envi- 
ronments. A large probe speed enforces a perfect Marko- 
vian behavior which matches the analytical curve (see OH)- 
Fig. EJb) shows the concurrence at a given time t a as a func- 
tion of the distance for three different entangled states: two 
Bell states and a cluster state (see figure caption). For Bell 
states the concurrence is equal to the absolute value of their 
only non-zero off-diagonal element in the density matrix, and 
therefore Fig. ^provides direct information about the individ- 
ual coherences. The figure clearly shows the influence of cor- 
related collisions: coherences poi.io (appearing in are 
robust against correlated noise, coherences poo.n (appearing 
in \4> + )) are especially fragile under correlated noise, and the 
remaining coherences decay in the same way under correlated 
or uncorrected noise (hence, the weak distance dependence of 
\G)). From Fig. |2b) it is also evident that the environments 
become more and more independent as d increases. In the 
limit of large distances (i.e. independent channels) the con- 
currence of the Bell states and the cluster state are related by 
C G = max{0, i(-l + 2 V /CS + C BeU )}. 

In this paper we studied the stochastic generation of entan- 
glement in a spin gas where the classical kinematics of the par- 
ticles drives the quantum state of the many-body system. This 
provides a novel scenario to study relationships between clas- 
sical thermodynamical variables and properties of the quan- 
tum state of the system such as quantum correlations, rate of 
entanglement generation, saturation values and clustering ef- 
fects. We have shown how to formalize such systems under 
quite general conditions and examined in some detail the cases 
of the Boltzmann and the lattice gas. We have fully charac- 
terized the bipartite aspects of the entangled TV-body state at 
long times, showing that maximum entanglement is asymp- 
totically reached for all possible bipartitions. The formal- 
ism has also enabled us to treat exactly a microscopic model 
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FIG. 2: (a) Concurrence of two probe particles prepared in a maxi- 
mally entangled Bell state = l/v / 2(|00) + |ll)) in a 100x3200 
lattice gas with N = 8 x 10 4 for g = O.877 and fixed probes 
(solid); and fast moving probes (dashed). Thin-solid curve cor- 
responds to analytical result obtained for the Markovian regime, 
(b) Concurrence of two probe particles at time t = 25rj _1 for 
different initial entangled states as a function of relative distance: 
\4> + ) = l/v^dOl) + 1 10)) (dashed-dark) , |^ + ) (solid-light), and 
\G) = l/\/2(|+0) + 1-1)) (dotted). 



for non-Markovian decoherence, or in other words, the time 
dependent maps describing the decoherence process of a set 
of probe particles that are in contact with the gas. We have 
seen that by changing the system parameters one can reach 
qualitatively different behaviors. The proposed models open 
the door to study entanglement dynamics in mesoscopic dis- 
ordered systems. The recent experimental breakthroughs on 
cold atoms in optical lattices I18i Il9l 12(11 provide promising 
perspectives for the experimental study of spin gases. 
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